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1 . INTRODUCTION 


An additional task of the human pilot of aircraft besides 
the primary task of being the controller of the aircraft is 
the monitoring task to promptly cope with troubles and sudden 
changes in the automated devices and the measurement 
instruments. It is quite important to become thoroughly 
familiar with the human function which detects such troubles 
and changes in the characteristics of the system and to 
express it using mathematical models. This is important for 
the designing of aircraft as well as other human - machine 
system. 

Research activities in the area of formation of models of 
human function in human - machine system are many. In these 
research activities, human is considered as servo mechanism. 
1-4) At the same time, research activities of human 
compensatory control that is seen at the time of sudden 
changes in the characteristics of the system are also carried 
out. 5) Yet, these research activities are not related to the 
formation of models of abnormal detection system of human 
being. Instead, they are thought to be the consideration of 
the model formation related to the control of the closed-loop 
system. 

Past research activities concerning human being from the 
standpoint of the detection of the system abnormality include 
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reports on experimental examination of human detection 
function. 6,7) On the subject of the model formation 
concerning the function of human detection, some of them 
express the human detection model of the monitoring task by 
the detection of gradual probability ratio applied to Carman * 
filter and integral calculus of the remaining difference of 
the filter. 8) The above model appears to be suitable in the 
case where abnormality appears as the bias change in the 
output. Yet, it (the model) seems to require further 
discussion if it is to be applied to the situation where the 
target system shows sudden changes or where sudden changes are 
made in the controlled element as human is carrying out the 
control within a closed-loop system. 

Authors conducted an experiment aimed at the secondary 
system on the subject of detection by human of parameter 
changes of dynamic system involving irregular noises as 
excitation input. 9) Based upon the result of the above 
experiment, authors have reported that the detection motion in 
the case where human is simply functioning as a monitor can be 
expressed using the most suitable assumption and the 
variance-ratio test. 10) 

In this research which deals with the detection motion at 
the time when the pilot control is in effect, above mentioned 
motion as the monitor is also included as the special case. 

In the same manner as the case of the monitor, the detection 
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time by the experiment (time between the change in the 
parameter and the time when the change is detected) can be 
connected to the signal being watched and dispersion of the 
time differentiated figure. . The mathematical model in this 
case consists of the control of the closed-loop system and the 
detection of abnormality. Model formation of the former is 
done by the optimal controller through the Weiner's method. 
Model formation of the latter is done according to the same 
method as that of the monitoring tasks. For the purpose of 
discussing the adequacy of the above model, a digital 
simulation of abnormality detection is conducted through the 
use of a model in order to compare the result with the result 
of the experiment by human. 


2. THE OBJECT SYSTEM 

Fig. 1 shows the control task used here. The following 
secondary system are used as the controlled element G(s): 


0(s)= -r~^ — ; 

S + 2 {„WnS + (On 

S 7 > — 9 

& : &rnt 

(On ■ SWiHfilSC 

K : 


( 1 ) 


S : parameter of rapuras * exchange 

r a : damped ratio 

iOn. : specific frequency 

K : Gain coefficient 






Fig. 1 Control task 


human 

controller Yp(s) noi** 

output 

y 


The way the abnormality appears caused by mechanical 
trouble and others tends to vary. Some are drastic. Examples 
of the above are rapid change in the number of the system or 
polarity. Others are less drastic. Examples are bias change 
in output, rapid change in the parameter K,^^, or ajn, and 
change in characteristics of the excitation input u. In this 
study, the question of detecting as soon as possible only the 
parameter^ ft of the controlled element G(s) or onlycOn. or the 
rapid changes which take place concurrently is handled as a 
human conducts compensatory control so as to make the error e 
zero under the influence of irregular noise u. Such sudden 
change is thought to be equivalent of a mechanical trouble 
developed in the safety devices of aircraft. Human 
transmission coefficient Yp(s) is explained in chapter 4.1. 

5. EXPERIMENT 

In conducting the experiment, controlled element G(s) was 
realized on an analogue computer. Tested personnel was 
provided with the error signal e as visual information. This 
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was shown on an oscilloscope as a horizontal line which moves 
toward up and down directions. The tested personnel was 
directed to maintain the horizontal line at the position as 
close to the standard line as possible. The operating device 
with a joy stick is structured to generate output voltage in 
proportion to the rotation angie of the joy stick. The figure 
of the control gain Kc was set at 0.1 volt/ deg or twice this 
figure. (The test result shown in this essay describes only 
the case of 0.1 volt/ deg.) Furthermore, the tested personnel 
was directed to press a push switch when he (she) made a 
judgment that some changes have occurred in the parameter 
through the observation of the error signal. Under the normal 
condition (shown in the letter n),^ rt = 0.7 and u>a = 4 rad/ sec 
at all times. Under this condition, only the parameter , 
only the parameter or the two parameters were changed. The 
excitation input u is the output of the noise generator. Its 
spectrum density is fixed up to 500 Hz. It is Gaussian signal 

» J 2 

whose figure is 7.24 x 10 V / Hz. In all experiments, the 
gain K of the controlled element is set at 1. 

The experiment was conducted on two subjects (two 
graduate students: 0 is 23 years of age, and H is 26 years of 

age.). Eight series of tests as shown in Tab. 1 were 
conducted on each subject. The experiment lasted 
approximately an hour and a half a day. The eight series of 
tests were completed within two days. In conducting each 


series of tests, the order of the test was from the ones with 
greater parameter changes to the ones with less parameter 
changes. The time when changes in the parameter occur is the 
optional time between 10 to 25 seconds after the commencement 
of each test. Data collection was initiated after the tested 
personnel became adequately accustomed to the parameter 
changes. Measured data are in the form of time (detection 
time td) between the parameter change in the controlled 
element and the time when the subject of the test presses the 
push switch. At each test point, approximately ten trial 
tests were conducted. If the push switch was not pressed 
within fifteen seconds of the parameter change, then such data 
were considered as "miss". If the push switch was pressed 
before the parameter change, then such data were classified as 
"false" and removed from the measurement data. The result of 
the experiment shown later is the average figure of the data 
collected according to the above method. Among the detected 
time, the longest and the shortest ones were removed before 
such data were averaged. 

Tab. 1 t'orin of changes in parameters 
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. Fig. 2 and Fig. 3 show the test result in which£n. and u)* 
make changes alone independently. Fig. 2 shows the 
measurement data of the detection time in which onlyjk changed 
(series 1), 2)). Fig. 3 shows the measurement data of the 
detection time in which only {*!«. changed (series 3), 4)). 

Circles with a dot show the test result, while blackened 
circles indicate the average detection time in the case of 
monitoring tasks. They show a little longer time than the 
average detection time in the control tasks. The mark r 
indicates the detection time as tne horizontal line on the 
oscilloscope was made to jump by 1.5 cm. It is considered as 
being equivalent of the reaction time of the test subjects. 

The total sum of the data used (including the "miss" and 
"false" data) is shown in TND. In the same figure, the solid 
line curves indicate the relation of the ratio of the 
detection time td and the average of the output y before and 
after the parameter change and square average figures of y 
using the following formulas (2) and (3). 



( 2 ) 



( 3 ) 


In the above formulas , tfy 1 anddy*indicate the square average 
figures of y and y. The letter alpha indicates the figure 
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after the parameter change. The figures of K^y and K^y can be 
obtained through the minimum square method from the test data. 
They are shown in the figure. (The ones within parentheses 
indicate the ones in reference to the monitoring tasks.) The 
curves in the figure shown in the formulas (2) or (3) are in 
reference to the controlled tasks. Yet, they are all in 
agreement with the test data, which indicates that the 
detection time by human can be well related with the ratio of 6y Z 
and<5y before and after the parameter change. 

Fig. 4 Detection time curve 



® 4 @ < 

3 > I- o — 7 1 — . |+ * £ 

Solid iines indicate the controller. Dotted lines 
indicate the monitor. Dots indicate the test points. Figures 
indicate the detection time (sec). Figures in parentheses 
show the case of the monitor. 
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Fig. 4 shows the result of the series 5) - 8) in which £n, 
and 60n are simultaneously changed. The horizontal axis 
indicates the ratio of<5y*before and after changes; vertical 
axis indicates the ratio otdy\ The dot marks indicate the 
test points. Test points of the series l) - 4) are also shown 
here. The numbers attached to the test point of the series 5) 
- 8) indicate the average detection time. (Figures in 
parentheses indicate the case of monitoring tasks.) The 
curves in the same figure indicate the average detection time 
using the following formula (4). 



K1 2 and K2 2 were determined through the minimum square 
method using the entire test data. The solid lines indicate 
the case of the control tasks, while the dotted lines indicate 
the case of the monitoring tasks. The direction toward which 
and become smaller and larger together, particularly the 

direction of the changes of the series 1) and 2) tests, 
compared with other directions is considered as the direction 
with longer detection time, in other words, the direction that 
is difficult to detect in spite of the fact that the change 
ratio of dyf and<Jy„ J is great. On the contrary, the direction 
of the series 5) and 8) is the direction with short detection 
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time in spite of the small change ratio. In other words, it 
is considered as the direction that is easy to detect. As the 
monitoring tasks and the control tasks are compared, the 
detection time is shorter in the case of the controller rather 
than the observer in the direction that is considered 
difficult to detect. In the direction that is thought to be 
easier to detect, no significant changes are recognizable in 
the detection time. 

Fig. 2-4 are in reference to the tested personnel H. 
Similar result has been obtained in the case of the tested 
personnel 0. For further detailed experimental data, refer to 
the reference literature 9). 

4. HUMAN DETECTION MODEL 

From the result of the experiment described in the 
previous chapter, it became apparent that human detection time 
was related to the ratio of the square average figure of y and 
y (equal to -e and -e here) before and after the parameter 
change. Accordingly, numan is thought to utilize the monitor 
volume as well as its differential figure at the time of 
structuring human detection model. The most suitable estimate 
is presumably provided for the purpose of obtaining the signal 
e and e which are used at the time of detecting the parameter 
change. At that time, human is thought to utilize his own 
output (operational volume m for the control) as known 


(established) signal. In the case of the monitoring, m 
becomes equal to zero. 

It is necessary to consider human control characteristics 
in the closed-loop system in order to calculate the output m. 
Human is assumed to conduct the control of the controlled 
element prior to the parameter change until he notices the 
change even if sudden changes occur in the controlled element. 
The above assumption is necessary in order to prepare the 
model of the control characteristics. As for said model, the 
optimal controller model through Weiner's * method is adopted. 
It is considered in the bibliography 4). 

As human detection model, the model as shown in Fig. 5 is 
proposed. The part of the most suitable assumption is modeled 
by Carman * filter, and the parameter change detection part is 
modeled through the use of the variance-ratio test. These 
parts, other than the influence of the human output m, are 
considered the same as the model formation 10) of the motion 
as the monitor. 

4.1 Optimal controller (within Fig. b) 

In this section, human transmission coefficient 
Yp(s) which would minimize a certain evaluation standard as 
the controlled element G(s) is provided is obtained. The 
problem in the bibliography 4) is solved as follows. The 
wasted time and the delay time of the nerve - muscle system 
that is expressed by 1/ (TNS + 1) among the things which show 
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human basic control motion is considered as the same in human 

who are matured up to a certain degree. This is not something 

that can be changed intentionally. In this section, wasted 

time as well as the delay between the nerve and muscle system 

-r$ 

are expressed by e together and left within the human 
transmission coefficient. Next, human control motion that is 
often called the proportional motion, the differential motion 
and the integral motion is the part which shows adaptability 
to the characteristics of the controlled element or the 
external disturbance. In other words, these are the parts 
that human can intentionally change. These parts are 
determined through the optimal controller theory using the 
following evaluation standard. 

_ i 5 .? 1 

In the above formula, the lines on top of letters e and m 
indicate the average figure. The second section on the right 
shows human labour. In general, m is considered as the item 
to show human labour. Yet, it is assumed that no additional 
labour is needed in order to maintain it at a fixed level 
after completing a fixed operation in the case where 
instruments without recovery power are used. 

Fig. 5 Human detection model 
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It is more proper to use m* which is equivalent of the 
average kinetic energy as the one to indicate human labour. 
Adequacy of the evaluation standard shown in the formula (5) 
is discussed after comparing it with the test result obtained 
by the National Aerospace Laboratory. 

According to the method described above, optimal 
human transmission coefficient Yp(s) is given according to the 
following formulas. 


Y ( s \— _ 1 F(s) 

Yv(s) G(s) 1 -«T"F(s) 






e r *G(s)G(-s) 


I 


Z(- s) 

G(s)G(-s)-*,V=Z(s)Z(-s) 


(6) i 

| 

(7) j 

( 8 ) 


Z(s) is a function structured by a pole of the ieft side half 
of the surface (L. H. P.) of the complex surface and the zero 
point. [•] + indicates the section which has a pole only at 

L. H . P. 

tig. b snows human transmission function Yp(s) 
obtained through the above described method in the form of 
board line chart when the controlled element G(s) is given by 
the formula (l). The result is equal to the result in 
reference to G(s) prior to the change in the parameter. As 





the figure of km 2- , kinetic movement in the vicinity of the 
roll axle by movable simulator is simulated in the test. This 
is a little different from the one described here. As an 
experimental attempt, km 1 = 0.05 sec 1 that is compared with 
the test result in the bibliography 4) was adopted. As the 
figure of r, the instance of 0.15 - 0.3 sec is shown here. 

Yet, based upon the following reason #1 and #2 , = 0.25 sec 

has been estimated in the digital simulation that follows. 
[Reason #1 ] The figure of T that is gained as the result of 
the experiment is 0.238 sec in the case of the test subject H 
and 0.306 sec in the case of the test subject 0. [Reason #2] 
The estimated figure of f in the bibliography 4) is almost 
0.25 sec. Incidentally, the indentation in the vicinity of 
the frequency 4 rad/ sec shown in the Gain's characteristics 
means the reduction in the pilot's gain in the vicinity of the 
specific frequency of G(s). 

4.2 Optimal estimator (in Fig. 5) 

Human internal observation noise v is determined as 
the Gaussian white noise. The excitation input u into the 
system uses the Gaussian signal. Also, it is assumed tnat 
human are adequately used to the signal prior to the parameter 
cnange and are well aware of these statistical 
characteristics. Considering tne fact that the opposing 
system is linear, the part of the optimal estimator of the 
detection model will become the Carman * filter assuming that 
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human are conducting the optimal estimation in the sense of 
minimizing the average square errors. The formulas for the 
optimal estimator by the Carman * filter are reported in 
various literature. Therefore, they are omitted here. 

fig. 6 A board line figure of human transmission function 



#6 0 Ara©£3SM#t<D#-KI*0 


4.3 Variance-ratio test (in Fig. 3) 

When the average figure of a base group N1 and N2 
(which are different from each other) is unknown, 
variance-ratio test to determine wnether some variance exist 
between^) anddj is considered here. The size of the sample 
that is extracted from each base group is determined at nl and 
n2. At this time, equitable estimate value ul and u2 in 
reference todi and d* are used to calculate the formula (9) 
shown below. 
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( 9 ) 



Then, L follows the F distribution of F (nl - 1, n2 - l) under 
the hypothesis ofO, = 0* . (nl - 1 and n2 - 1 indicate the 
degree of freedom.) F( alpha : m, n) is defined as the figure Atf- 
for the probability variable Z following the F (m, n) 
distribution to satisfy P(ZiAn) = alpha. Then, if either one 
of the formulas (10) or (ll) shown below is materialized, then 
the hypothesi s 6* = d* is discarded. If none of the two 
formulas is materialized, then the hypothesis is not 
discarded. This is the variance-ratio test in mathematical 
statistics. Alpha becomes the danger ratio to discard the 
hypothesis at the time when the hypothesis is correct. 

Fig. 7 shows the method to apply the above method to time 
series signal. The two parts among said signal that will 
become the target of the variance-ratio test are Nl and N2. 

In other words, it is estimated that human conduct 
variance-ratio test by extracting nl ana n2 each samples from 
Nl and N2 part. In order to detect the parameter change, 
human observes this signal. The actual point of the 
observation is at the right edge of N2 . Considering the fact 
that human are quite familiar with the statistical 
characteristics under the normal condition, it may be ail 
right to think that Nl is positioned at the beginning of the 
signal. It was thought that only the part of N2 was moved to 
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(10) 


L2F(a/2: »,-l, n,-l)&F„, t 

the right along with the passage of time. *, 

HL^F(a/2 : m-1, (n) 

5. DIGITAL SIMULATION BY DETECTION MODEL 

Digital simulation was conducted using human model 
proposed in the chapter 4. The condition of the set up was 
made identical to that of the experiment. The sampling period 
was set up at 0.05 second. It was so decided considering the 
fact that the reading accuracy at the time of the measurement 
of the detection time during the experiment was 0.05 second 
and thatu>«= 4 rad/ set (accordingly, the specific period is 
1.57 sec). As for the human internal observation noise, it 
was determined so that the ratio of the variance of e(^e*)and 
the variance of v(<fv)was equal tode Z /df - 36. 8) The length of 
the signal extracted from the part of N1 can be anything as 
long as it is long enough to determine the statistical 
characteristics. Cons idering Wn = 4 rad/ sec, it was 
determined as follows: N1 = 10 seconds. As for the moving N2 

part, it can also be determined according to the same manner 
as Nl. Yet, if this part is too long, then the detection 
becomes somewhat delayed. If it is too short, then it could 
cause some error because of its inability to fully grasp the 
statistical characteristics. Here, N2 is made equai to 4.5 
seconds. The time when the parameter becomes changed is 12.5 
seconds after the beginning of the test. 

Fig. 8 to Fig. 13 show the result of the simulation of 
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the detection time through the use of models. Fig. 8 and 9 
show the change of solely^ , while Fig. 10 and 11 show the 
change of solely Ct)* . Fig. 12 and 13 show the result in which 
both£ ft and change simultaneously (series 6), 7) ). In 
these figures, detection time obtained by the experiment using 
human is also shown using A mark. Fig. 8, 10 and 12 show the 
detection time obtained from the variance-ratio of the optimal 
estimate figure e' of e. Fig. 9, 11 and 13 show the detection 
time obtained from the variance-ratio of the optimal estimate 
figure t of e. The sample numbers nl and n2 to be extracted 
from the part of Nl and N2 only need to be sufficient to know 
the statistical characteristics of tne signal. It should be 
determined based upon the frequency band of the signal. The 
number of the sample usually affect the threshold figure (for 
detection) as well as detection time or generation of "false". 
Here, the result in the case of (nl, n2) = (200, 90), (100, 
45), (66, 30) is shown. It is the result obtained when alpha 
= 0.01. The threshold figure F crit for the detection at that 
time is shown in Fig. 2. Each result is the average of the 
result of 100 separate simulations whicn are correctly judged. 
Comparing the above result with tne result of the experiment 
marked with.A, it is obvious that they tend to snow decent 
agreement. In the case of the change of alone, the 
detection time in the case where ul« 2 changes toward the higher 
(.larger) direction is shorter than the detection time in the 


case where the same changes toward the lower (smaller) 
direction. This is explained according to the following 
manner, in the domain where is high (large), as it has been 
reported previously that a pilot will begin the control after 
noticing <»)n3) , <^«*can be more easily detected. This tendency 
also appears in the detection time obtained from e (fig. 11). 
It is safe to assume that human is paying more attention to e. 

In the case of the series 5) and 8), it was difficult to 
make the detection through the simulation according to this 
model. This is the case in which both and became larger 
or smaller. Human easily detected the change in the parameter 
in this case in spite of the fact that the change ratio of the 
variance of e and e was rather small. The above fact 
indicates the possibility that human commonly use the key 
information other than the variance-ratio in order to detect 
the changes in parameter. However, as far as other direction 
is concerned, human perception is directed only toward the 
change ratio of the variance. Thus, fine result is obtained 
as shown here. 

The influence of the control is not significant as the 
result of the simulation is compared witn the result of tne 
monitoring tasks in tne bibliography 10). It is particularly 
so when all the parameters are regarded as being the same 
(equal). Assuming that the number of samples, nl and n2, from 
the Nl and N2 parts is greater in the case of the controlled 
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tasks than in the case of the monitoring tasks (for example, 
monitor (nl, n2) = (100, 45) and the controller = (200, y0) ), 
approximately the same degree of difference as the difference 
of the experimental figure of the detection time in the both 
cases will appear. This can be interpreted that human tend 
to be more attentive in the case of being a controller rather 
than being a monitor. In this study, the data length of N2 
part is calculated as being fixed. It is assumed that the 
detection characteristics may be changed by changing it. In 
this essay, no discussion has been made concerning the method 
on how to select these parameters. Consideration of the above 
discussion may be brought out in the future. 

6. CONCLUSION 

In this study, an experimental investigation of detection 
actions of airplane pilots was conducted in the case where the 
system abnormality can be detected as the variance changes of 
the output, based on the above study, a detection model was 
proposed. It became clear that the pilot's detection time 
obtained through experiments could be connected with tne 
variance-ratio of the signal being watcned and its time 
differentia] figure. Based upon the result of the 
experiments, the pilot's detection model is expressed in the 
form of the two parts - the optimal control for the 
closed-loop system control and the detection of the parameter 
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changes by Carman * filter and the variance-ratio test. 

Through the comparison of the result of the digital simulation 
by models and the test result, it is safe to state that the 
model being proposed here adequately and qualitatively 
expresses pilot's detection motions as a controller or a 
monitor with the exception of the special case such as the 
series 5 ) and 8) . 

The questions yet to be solved may include tne mutual 
relation between e and e which are utilized by pilots who are 
attempting to detect changes. In other words, the future 
consideration must include the question on which one must be 
used, or how much of it should be used, or now it should be 
used... etc. etc. In addition to that, selection method of the 
parameters in the use of the models must be discussed through 
further tests and experiments including the research 
activities on the subject of the difference between the case 
of the monitoring tasks and the control tasks. It seems that 
consideration of the control standard other tnan the 
variance-ratio is also important. 

[Translator's note] 

Words marked with * indicate phonetically translated 
Japanese words. 

Bibliography 

l) McRuer, D.T. and Jex, H.R.: A Review of Quasilinear Pilot 

23 

... : - • 4 .... . ■ 


Models, IEEE Trans. Hum. Factors Electron, HFE-8 (1967), pp. 
231-249. 

2) Kleinman, D.L. and Baron, S.: Manned Vehicle Systems 

Analysis by Means of Modern Control Theory, NASA CR-1753 June 
(1971). 

3) Kimura, H., Goto, N. , Morizumi, N. , Fukuda, S. and 

Kitagawa, G. : On the Manual Control of a Second-Order System, 

Theor. Appl. Mech. , 25 (1977), pp. 3-11. 

4) N. Sumiyoshi , J. Kumamoto, S. Goto: Consideration on 
optimal kinetic characteristics of pilots, Ningen kogaku, 25 
(1979) , pp. 87-95 

5) Young, L. R.: On Adaptive Manual Control, Ergonomics, 12 
(1969), pp. 635-674. 

6) Curry, R. E. and Govindaraj, T.: The Human as a Detector 

of Changes in Variance and Bandwidth, Proc. 13th Annual 
Conference on Manual Control (1977), pp. 217-221. 

7) Curry, R. E. and Govindaraj, T. : The Psychophysics of 

Random Processes, Proc. 12th Annual Conference on Manual 
Control (1976), pp. 50-59. 

8) Gai, E. G. and Curry, R. E. : A Model on the Human Observer 
in Failure Detection Tasks, IEEE Trans, on System, Man, and 
Cybernetics, SMC-6, No. 2 (197b), pp. 65-94. 

9) N. Morisumi , S. Goto, H. Kimura: Detection of system 
changes by pilots: Experiment, Kyudai Kogaku Shuho, 56 (1983), 
pp. 71-77 


24 


10) N. Morisumi, S. Goto, H. Kimura: Characteristics of system 
change detection by human and its models. System and control, 
27 (1983), pp. 532-537 


Tafc Z 


* 2 * 



F«rll^ = 

0.01 

i 

! 

m 

n 

Fern 

199 

89 

1.621 

99 

44 

2.002 

65 

29 

2.432 

89 

199 

1.564 

44 

99 

1.882 

29 

65 

2.165 


F-.S. 7 





3 • 


25 






